Gauge-invariant critical exponents for the Ginzburg-Landau model 
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The critical behavior of the Ginzburg-Landau model is described in a manifestly gauge-invariant manner. The 
gauge-invariant correlation-function exponent is computed to first order in the 4 — d and 1 ^-expansion, and 
found to agree with the ordinary exponent obtained in the covariant gauge, with the parameter a = 1 — d in the 
gauge-fixing term (d^A^) 2 /2a. 
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Despite being one of the most studied field-theoretic 
models in theoretical physics, the critical behavior of the 
Ginzburg-Landau model is still poorly understood due to non- 
trivial gauge properties. The model is defined by the Hamil- 
tonian 



(1) 

where F^ v = d^A v — d v A^, with \i, v — 1, - ■ -d, e and 
m are electric charge and mass, and A parametrizes the self- 
interaction. The last term with parameter a fixes the gauge. 
The phase transition occurs where m 2 changes sign. The com- 
plex field (j> plays the role of an order field which has a nonzero 
expectation value in the ordered phase. 

Apart from the standard field-theoretic interpretation, the 
Ginzburg-Landau model can be equivalently understood as 
describing a random tangle of intertwined electric current 
loops of arbitrary length and shapeoB. In the normal state, 
only a few current loops are present due to a finite line tension 
0. At the critical temperature T c , the tension vanishes and the 
current loops become infinitely long. An important character- 
istic of these geometrical objects is their fractal dimension D, 
which at T c is related to Fisher's critical exponent 77, determin- 
ing the anomalous dimension of the order field (see below). In 
the absence of gauge fields, this exponent manifests itself in 
the power behavior of the correlation function 



G(x-x') = (<p(x)f(x')) 



(2) 



at the critical point as being G(x) ~ 1/x 2+n *. The free 
theory has G(x) ~ l/x d ~ 2 , corresponding to 770 — 0. A 
nonzero value of the critical exponent 77^ implies that the di- 
mension of cf> deviates from the canonical, or engineering di- 
mension (d — 2)/2. 

For a particular value of the Ginzburg-Landau parameter 
kgl, defined by the ratio Kq L = e 2 /A, the Hamiltonian ([[]) 
is also a dual description of an ensemble of fluctuating vortex 
lines of arbitrary length and shape in superfluid heliumlil, in 
which case is a disorder field. 

In this note, we wish to clarify the properties of this impor- 
tant exponent, which has been controversially discussed in the 
past, and very recently also in the context of quantum electro- 
dynamics (QED) — the fermionic counterpart of the Ginzburg- 
Landau model (see Refs. and references therein). 
The poor understanding of this exponent is because in a gauge 
theory, the correlation function (||) depends on the gauge pa- 
rameter a in Eq. (m. This has led to a severe theoretical puz- 



zle. In ordinary local quantum field theories without gauge 
fields, one can prove that it must be greater or equal to zero. 
In contrast, renormalization group studies have always pro- 
duced negative (albeit gauge-dependent) values, starting with 
the historic paper by Halperin, Lubensky, and Mail. In gauge 
theories, the proof of a non-negative 77 is not applicable due 
to the nonlocal nature of the gauge-invariant correlation func- 
tion, as has been understood only recently!]. 

We first consider the model close to the upper critical di- 
mension in d = 4 — e dimensions, and extend the Hamiltonian 
([]]) for later discussions to contain n/2 complex fields with an 
0(7i/2)-symmetric self-interaction. 

A one-loop perturbative treatment of the Hamiltionian ([[]) 
yields the first term in the e-expansion of the critical expo- 
nents. For the exponent 77^, the well-known result iaZlH: 
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where e 2 = 487r 2 e/n is the value of the charge at the fixed 
point. For an infrared-stable fixed point to exist in the two- 
dimensional space spanned by e and A, n must satisfy n > 
12(15 + 4\/l5) « 365.9 to first order in the the e-expansion. 
These results are for the massless model (m = 0). To avoid 
infrared divergences, Feynman diagrams are evaluated at fi- 
nite external momentum k. Being the only scale available, k 
is used to remove the dimension from dimensionful parame- 
ters, for example, e 2 = e 2 n d ~ 4 . Most calculations reported in 
the literature are performed in the Landau gauge (a = 0) for 
which0770 — > — 18e/n. 

Since the correlation function (||) is not gauge-invariant, it 
is not a physical quantity. It is therefore not surprising to find 
that the critical exponent 77^ depends on the gauge parameter 
a. As long as the gauge is fixed by the last term in Eq. ([l]), 
the correlation function is nevertheless well-defined and, as, 
following Ref. JTo|], can be verified explicitly to first order in 
e, the critical exponents satisfy the usual scaling laws: (3^ — 
\v{d - 2 + 770), dv = 04,(64 + 1), with 77^ given in Eq. (||). 
Since the correlation length exponent v is gauge-independent, 
ft,/, and 64, depend on a. The exponents 0$ and 5$ specify the 
vanishing of the order parameter ((f)), respectively for T — > T c 
from below and when an external field coupled linearly to <p 
tends to zero. 

A physical correlation function must be gauge-invariant, 
i.e., invariant under the combined transformations <j>(x) — > 
exp[ieA(x)](j>(x), A fl (x) — > A fl (x) — d fl A(x). Such a cor- 
relation function is obtained by including a path-dependent 
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Schwinger phase factor in Eq. ([|), forming 



G(x-x') = ^(x)^(x') 



(4) 



where the line integral extends from a; to x', and the av- 
erage, denoted by angle brackets, is taken with respect to 
the Hamiltonian (jl|). At the critical point, we now expect a 
power behavior G(x) ~ l/x d ~ 2+mi , where in contrast to 
?70, the exponent t/qi has no a-dependence. The exponential 
in Eq. (Q) can alternatively be written in terms of an exter- 
nal electric current line J^{z) as e _1 ^ d where 

x' 

Jfj,{z) = e f dx^5(z — x) is a delta function along the 
path from x to sc' which satisfies the current conservation law 
dftJ^z) = e8{z — x) — e5(z — x'), with a current source at 
x and a sink at x' . In Refs. [ p"lp2[ ], the gauge-invariant cor- 
relation function (Q), with its external current line, was stud- 
ied in d = 3 and found to behave differently in the normal 
and superconductive state. In the normal state, where the line 
tension 9 of the current line was shown to be finite, this cor- 
relation function decreases exponentially for large separation, 
G(r) ~ e~ 6r , with ry, = x'^ — x^ being the distance vector. In 
the superconductive state, on the other hand, the line tension 
vanishes and the correlation function was found to behave in- 
stead as G(r) ~ exp(e 2 AL/47rr), with Al being the London 
penetration depth. Rather than tending to zero, the correlation 
function now reaches a finite value for large separation. The 
finite expectation value at infinite separation signals that the 
current lines have lost their line tension and have become in- 
finitely long. In the correlation function (Q) it manifests itself 
in an independence on the path over which the line integral is 
taken. Only the endpoints of the line connecting x and x' are 
physical. The exponent in the correlation function contains a 
Coulomb-like interaction between these endpoints. (Note that 
the combination is independent of the electric charge.) 

To compute the exponent t]gi to first order in the e- 
expansion, the gauge-invariant correlation function ([I]) is ex- 
panded to order e 2 . Then, using Wick's theorem, three conti- 
butions are obtained besides the lowest order 



G(r)=G + T +T 1 +T 2 , 



(5) 



containing respectively no, one and two Schwinger phase fac- 
tors. The first contribution Tq is given by the integral 



d d zd d z 



G(x - z)d Zfl G{z'-x')d z >G{z'-z) D^{z-z') 



where the right minus left derivatives d Zll = d Zfl — dz^ oper- 
ate only within the square brackets, and is the correlation 
function of the gauge field in ([I]), with the Fourier compo- 
nents 



q A 



Sfiv - (1 - a) 
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In mometum space this yields 

2 f d d k d d q c ife -'' (?+§*!)„ 
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D^q-k/2). 



Since the fixed-point value of e 2 is of order e [see below 
Eq. (^J)], the integrals can either be evaluated directly in d = 4 
or using dimensional regularization and taking the limit d — > 4 
at the end. Either way gives with logarithmic accuracy 
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Adding the free scalar correlation function G(r), we obtain 
,a— 3 



G(r)+T (r)=G(r) 
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G(r)r- £ ^^ , 
(7) 



which reproduces the old result (£ 
The last term in Eq. (61), 
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G{r) d d zd d z'J fl (z)D fi „{z-z')J„(z'), (8) 



T 2 (r) 



factorizes from the start in a scalar and gauge part, with the 
second factor — which plays a central role in the study of 
the gaiig&|invariant order parameter of the Ginzburg-Landau 
modeOEl — denoting the Biot-Savart interaction between two 
segments of the external current line. Since the integrals in 
Eq. (|^) are line integrals, their value depends on the path cho- 
sen. We choose as integration path the shortest path connect- 
ing the two endpoints, i.e., straight lines and write 



T 2 (r) = — =-G(r) / dudu'r^D^v! - u)r]r u 



(9) 



after the reparametrization x h 



ur^, x^ = x^ + u r M , 
- x v and < u, vl < 1. 



with a fixed distance vector r u = x,, 
The integrals are easily evaluated following Ref. |]3|], with the 
result 



T(d/2-l) 



(3-d)(4-d) 47r d / 2 
which for d near 4 yields 
,3- 



l-i(l- a )(3-d) 
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G(r)r^, 
(10) 

(11) 



Due to the appearance of the logarithm multiplying G(r), this 
gives the contribution 



(12) 



to the Fisher exponent. When both contributions obtained 
so far are subtracted rather than added, one obtains a result 
(which happens to be zerol independent of a. As first noted 
in the the context of QEDd, this is because the combination 
V4> ~ V2 characterizes the correlation function 



(13) 



(<p(x)$(x')){exp\ie f dx^A^x) 
which is gauge invariant. 
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Next, the third, or mixed term in Eq. 
?ral 



b, given by the inte- 



e / d d zd d z' 



(14) 

is evalutated. We expect an a-dependent contribution that 
precisely cancels the dependence on the gauge parameter 
found in Eqs. (|3|) and (|l2|). To extract the term of the form 
G(r) ln(r), we use the approximation, cf. Ref. [Q], 

Ti(r) w eG(r) f d d z d d z' [d Zfi G(z-x')-d Zfi G(x-z)] 



(15) 



x J„(z')Du V (z-z'), 



valid with logarithmic accuracy. Both terms in the square 
brackets give the same contribution. Partially integrating 
this expression and using the identity (x^x ll /x 4 ') = 
[(3 - d)/2]d v x~ 2 in d = 4, we obtain 



a 



Tx(r) = ~e 2 ^G(r) / d d zG(z-x') 



(z - x) 2 



(16) 



giving (e 2 a/ Air 2 )G{r) ln(«;r) and thus a contribution to rj 



a 
4^2' 



(17) 
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As expected, this contribution precisely cancels the in- 
dependence in Eqs. (fy and (|l2|). More specifically, we obtain 
for the manifesly gauge-invariant correlation function 

~2 3 36 
Vgi = V<t> + Vi + V2 = -e — — = e. (18) 

This value for tjqi is twice that for 77^ obtained in the Landau 
gauge (a = 0). Both results coincide, however, when a — 
-3. 

In the current loop description, the critical exponent ijq^ 
termines the fractal dimension D of the current lines via 
D = 2 — r/Qi- With t?gi < 0, the fractal dimension is 
larger than that of Brownian random walks for which D = 2, 
implying that the current lines are self-seeking rather than 
self-avoiding, which makes them more crumpled than Brow- 
nian random walks. Although higher-order corrections may 
well change the sign of ?7gi> nothing in the cpptext of the 
Ginzburg-Landau model forbids negative valuesELj, provided 
?7gi > 2 — d, or D < d. In the limiting case D = d, the 
current lines would be completely crumpled and fill out all of 
space. 

Instead of an e-expansion, we may compute the gauge- 
invariant critical exponent 77GI nonperturbatively in the limit 
of a large number n of field components. Then rj can be ex- 
panded in powers of 1/n for all 2 < d < 4. 

The leading contribution in 1 jn generated by fluctuations in 
the gauge field is obtained by dressing its correlation function 
with arbitrary many bubbkdnsertions, and summing the entire 
set of Feynman diagramsLa The resulting series is a simple 
geometrical one, which leads to the following change in the 
denominator of the pref actor in the correlation function (0): 



e — 



2 (d- 1) 



(19) 



where the second term dominates the first one for small q in 
2 < d < 4. In Eq. (Q, c(d) stands for the 1-loop integral 



G(x-z)d^G(z-x') Uz^D^z-z 1 ), c(d ) = 



1 



(27r) d k 2 (k+p) 2 



T(2~d/2)T 2 (d/2- 1) 
(47r) d / 2 r(d - 2) 



(20) 

where analytic regularization is used to handle the ultraviolet 
divergences—Tja leading order in 1/n the value of r]^ for 2 < 
d < 4 readsBO 



2 4-d - (d- l)[4(d- 1) 

v<t> = -- 

n 



(4w) d / 2 c(d)r(d/2 + 1) 



(21) 



which depends on the gauge parameter a. For d = 4 — e, 
this result reduces to Eq. (||) obtained to first order in the e- 
expansion. The gauge dependence of 77^ is not always obvi- 
ous from the results quoted in the literature as often a spe- 
cific gauge is chosen from the start, for example, the Landau 
gaugeQ, where — * — 40/7r 2 n for d = 3. 

The term (Q) with the modified gauge-field correlation func- 
tion can be evaluated as before. To extract the dependence on 
ln(r) it proves useful to replace q d ~ 2 with q d ~ 2 + s in Eq. (|T^), 
so that the gauge-field correlation function in the large-n limit 
becomes 



DM - 
and to let 8 

T 2 (r) = ~ 



2 d- 1 



1 



c{d) q a 



-2+8 



8nu - (1 - a 



(22) 



at the end. This leads to 

8(d- 1) 



n (47r) d / 2 c(d)r(d/2 - 1) 

with r 5 /8 — * 1/8 + m(r), and 

4 (d- l){d- 1 



1 



1 — a 
d-2 



a) 



n {4n) d / 2 c{d)T(d/2) 



(23) 



(24) 



As a check note that, when subtracted from ry^ given in 
Eq. (pi]), this yields an a-independent result characterizing 
the gauge-invariant correlation function ([jj]), 
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n (4ir) d / 2 c(d)T(d/2 - 2) ' 



(25) 



This expression is negative for all 2 < d < 4. Specifically, 

V4>-V2 = -8/n 2 n for d = 3 and -4e 2 /" + d{e 3 ) for d = 
4- e. 

To calculate the mixed term (|l4]), the gauge-field corre- 
lation is needed in coordinate space. Fourier transforming 
Eq. (22), we arrive at 



D^(x) 



2 d- 1 



1 



ne 2 c(d) (47r) d / 2 r(d/2) .x 2 -- 5 



3 + + (1 - a) 



(26) 



Proceeding as before, we obtain 
T l( r) = -2 16 «-V 



n (4ir) d / 2 c(d)T(d/2 - 1) 
1 



G(r) 



(27) 



d d zG(z 



(z - x) 



2-8 



-rjiGir) ln(r) 
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with 



a 8(d-l) 



n (4ir) d / 2 c(d)T(d/2) ' 



(28) 



which, being proportional to a, should cancel the a- 
dependence in -q<p and 772. And indeed, as grand total we find 
a result , , ,2 



VGi = rj<t>+ m+ f i2 = — 



(d 2 + 2d-6)r(d-2) 
nT(2-d/2)T 2 {d/2-l)T{d/2) 



(29) 



which is independent of the gauge parameter a. Remark- 
ably, for the d-dependent gauge choice (of which our a = — 3 
found in the e-expansion is a special case) the 77^ of Eq. ( pT| ) 
coincides with 77G1 to this order in 1/n. With this gauge 
choice, the trace of the gauge-field correlation function van- 
ishes, _D MM (g) = 0. Since this observation does not depend on 
the matter part of the theory, we expect it to hold also in QED. 
And indeed, the valuescH for the the two 77 exponents obtained 
in first order in 1/n in d = 3 and also in the e-expansion co- 
incide when a — 1 — d. Although higher-order corrections 
might change this simple relation between the two exponents, 



we speculate that in first order the gauge choice a = 1 — d 
provides a shortcut for obtaining the gauge-invariant result of 
other quantities such as the effective potential and mass renor- 
malization. 

The expression (^9|) is negative for all 2 < d < 4, with 
VGI = -72/n 2 n for d = 3 and -36e/n + 0{e 2 ) for 
d = 4 — e. The latter result is in accord with Eq. ( |l8| ) ob- 
tained in perturbation theory. We repeat that in the context of 
the Ginzburg-Landau model, negative values are allowed, pro- 
vided ?7qi > 2 — d, so that the fractal dimension D = 2 — t/qi 
of the current lines is smaller than the dimension of the em- 
bedding space. A negative value merely indicates that the cur- 
rent lines are self-seeking rather than self-avoiding. 
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